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Positive-strand RNA virus genomes are translated into polyproteins that are processed by viral proteases to yield functional in- 
termediate and mature proteins. Coronaviruses (CoVs) carry genes that encode an nsp5 protease (also known as 3CLpro or 
Mpro) responsible for 11 maturation cleavages. The nsp5 structure contains two chymotrypsin-like domains (D1 and D2) anda 
unique domain (D3), and forms functional dimers. However, little is known of interactions or communication across the struc- 
ture of the protease during nsp5 activity. Using reverse genetic mutagenesis of the CoV murine hepatitis virus (MHV) nsp5, we 
identified a new temperature-sensitive (ts) mutation in D2 of nsp5 (Ser133Ala) and confirmed a ts residue in D3 (Phe219Leu). 
Both D2-tsS133A and D3-tsF219L were impaired for viral replication and nsp5-mediated polyprotein processing at the nonper- 
missive temperature. Passage of tsS133A and tsF219L at the nonpermissive temperature resulted in emergence of multiple sec- 
ond-site suppressor mutations, singly and in combinations. Among the second-site mutations, a D2 His134Tyr change sup- 
pressed the ts phenotype of D2-tsS133A and D3-tsF219L, as well as the previously reported D2-tsV148A. Analysis of multiple 
CoV nsp5 structures, and alignment of nonredundant nsp5 primary sequences, demonstrated that ts and suppressor residues are 
not conserved across CoVs and are physically distant (>10 A) from each other, from catalytic and substrate-binding residues, 
and from the nsp5 dimer interface. These findings demonstrate that long-distance communication pathways between multiple 
residues and domains of nsp5 play a significant role in nsp5 activity and viral replication, suggesting possible novel targets for 


non-active site inhibitors of nsp5. 


Petar RNA viruses are responsible for prevalent and 
epidemic diseases in a wide range of vertebrate hosts, as well as 
new and emerging viruses, such as severe acute respiratory syn- 
drome coronavirus (SARS-CoV), West Nile virus, and Chikungu- 
nya virus. The rapid evolution, host species movement, and dis- 
eases of positive-strand RNA viruses demonstrate the need to 
develop novel strategies to prevent and treat present and new dis- 
eases caused by these viruses. A key determinant of positive-strand 
RNA viruses is the requirement for processing of translated poly- 
proteins by virus gene-encoded proteases. RNA virus proteases 
therefore have been high-profile targets for development of anti- 
viral agents, with most protease inhibitors targeted to active sites 
or substrate-binding sites (21, 26, 32, 35). However, due to the 
potential for viral escape mutants, it is critical to identify addi- 
tional noncatalytic, non-substrate-binding determinants of pro- 
tease activity as potential targets for inhibition that are less prone 
to development of resistance. 

To date, five CoVs have been shown to be associated with hu- 
man respiratory diseases of different degrees of severity: human 
coronavirus HCoV-229E, HCoV-OC43, HCoV-NL63, HCoV- 
HKU1, and SARS-CoV (10, 15, 30, 31, 44, 45). CoVs contain the 
largest known positive-strand RNA genomes, ranging from 26 to 
32 kb in length. Murine hepatitis virus (MHV) strain A59 is an 
established model for study of CoV replication and pathogenesis. 
The 32-kb genome of MHV contains seven genes, with the repli- 
case gene (22 kb) encoding 16 nonstructural proteins (nsp1 to 
nsp16) (Fig. 1A) (20, 25). The replicase gene is translated into 
polyprotein la (ppla; nsp1 to nsp11) or, via a ribosomal frame- 
shift, pplab (nsp1 to nsp16) (9, 25, 33). MHV encodes two papa- 
in-like proteases (PLP1 and PLP2) responsible for cleavages of 
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nsp1 to nsp3, and an nsp5 protease, also known as 3CLpro or 
Mpro, that mediates maturation cleavages of nsp4 to nsp16 and is 
required for virus replication (33, 50). 

The CoV nsp5 is a cysteine protease present in all known CoVs 
and is structurally similar to the nsp4 protease of distantly related 
arteriviruses (6, 33, 49). The crystal structure of nsp5 has been 
solved for divergent CoVs from every genus, including SARS- 
CoV, infectious bronchitis virus (IBV), human HCoV-HKUI1, 
and human HCoV-229E. Comparison of solved nsp5 structures 
demonstrates conservation of tertiary structure despite numerous 
differences in primary sequences (1, 2, 5, 46, 47, 49). The X-ray 
crystal structure of MHV nsp5 has yet to be determined; however, 
the structure of the closely related HCoV-HKU1 nsp5 (84% se- 
quence identity) has been resolved to 2.5 A (Fig. 1B and C) (49). 
The nsp5 proteases of all CoVs exhibit a three-domain structure, 
with domains 1 and 2 forming a chymotrypsin-like fold contain- 
ing the His41-Cys145 catalytic dyad and substrate-binding sites 
(Fig. 1B) (1, 2, 5, 28, 46, 47). In contrast, domain 3 is unique to the 
CoV nsp5 protease among chymotrypsin-like enzymes and also 
shows more divergence in both sequence and structural organiza- 
tion between CoVs. In vitro studies indicate that domain 3 is im- 
portant for stabilization of the chymotrypsin-like fold and may 
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FIG 1 MHV protease nsp5 structure and sequence alignment. (A) Genome and polyprotein processing of MHV. The replicase gene (gray) consists of two open 
reading frames (ORFla and ORFlab) overlapping at a ribosomal frameshift and is translated to yield two polyprotein products, ppla and pplab, encoding 
nonstructural proteins (nsps) nsp1 to nsp16. Maturation cleavages (arrows) are mediated by three viral proteases, papain-like proteases PLP1 (light gray) and 
PLP2 (dark gray) and nsp5 (black). M, matrix; E, envelope; N, nucleocapsid; RdRp, RNA-dependent RNA polymerase; Hel, helicase; ExoN, exoribonuclease; 
EndoU, endoribonuclease; O-MT, O-methyltransferase. (B) Modeled structure of a MHV nsp5 monomer based on the crystal structure of HKU1 (49) shows 
domain 1 (D1), domain 2 (D2), and domain 3 (D3). ts and second-site suppressor residues described by Sparks et al. are also identified (40). N and C denote the 
amino and carboxy termini. The dimerization interface is labeled, and functional residues discussed are shown in black and are labeled. (C) Sequence alignment 
of MHV-A59 and HCoV HKU] with nonconserved residues highlighted, catalytic residues boxed, and domain separations identified. 


also be important in mediating dimerization between nsp5 mono- 
mers (1, 28, 37, 38). Structural and biochemical studies demon- 
strate that nsp5 dimerization is required for proteolytic activity in 
vitro (11, 12, 37). 

Recently, amino acid residues that may regulate nsp5 activity 
but that are distinct from the active site cavity, substrate-binding 
pocket, or dimerization interface have been identified. Alanine 
substitution at SARS-CoV nsp5 Ser147, a conserved serine residue 
in MHV and HKU1, disrupts dimerization and impairs nsp5 pro- 
teolytic activity, despite being greater than 9 A from the dimeriza- 
tion interface (7). Our laboratory identified a temperature-sensi- 
tive (ts) mutation (tsV148A) in MHV nsp5 that impairs virus 
growth and nsp5 activity at 40°C (Fig. 1B) (40). Growth of the 
tsV148A mutant at 40°C resulted in emergence of second-site mu- 
tations (S133N and H134Y), which suppressed the tsV148A phe- 
notype. While V148 is adjacent to the catalytic C145, neither 
S133N nor H134Y have predicted direct interactions with catalytic 
or substrate-binding residues or identified pathways for propaga- 
tion of structural changes. Sawicki et al. identified a putative ts 
allele in nsp5 (F219L) at the base of unique domain 3, and distant 
from the dimerization interface, active site cavity, and substrate- 
binding regions (36). Finally, other studies have shown that nsp5 
activity may be altered by changes in replicase nonstructural pro- 
teins nsp3 and nsp10 (14, 41). Together, these findings support 
the hypothesis that residues distant from catalytic and substrate- 
binding sites are important for regulating nsp5 protease activity. 
However, mechanisms of communication between residues and 
regulation of protease activity are unknown. 

In this study, we tested the roles of MHV nsp5 domain 2 resi- 
dues $133 and H134, and domain 3 F219, on nsp5 activity. The 
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experiments identified a previously unknown MHV nsp5 ts mu- 
tation in domain 2 (tsS133A) and confirmed the ts phenotype of 
F219L, both of which result in profound growth and protein- 
processing defects at nonpermissive temperatures. Under passage 
at nonpermissive temperatures, a series of second-site mutations 
emerged that were able to suppress the ts phenotypes of tsS133A 
and tsF219L while being physically distant from the ts allele, the 
nsp5 active site cavity, and the dimerization interface. A single 
nonsynonymous mutation resulting in a H134Y substitution sup- 
pressed the ts phenotype of all three independent ts alleles in do- 
mains 2 and 3. The ts mutations reduced but did not abolish nsp5 
protease activity during virus infection following a shift to the 
nonpermissive temperature, while individual and combined sup- 
pressor mutations restored nsp5 activity to an extent that directly 
correlated with increased replication. These results demonstrate 
the presence of multiple interconnected long-distance communi- 
cation nodes in nsp5 and suggest novel mechanisms of regulation 
of nsp5 activity during CoV replication. 


MATERIALS AND METHODS 


Viruses, cells, and antisera. Recombinant wild-type (WT) MHV strain 
A59 (GenBank accession no. AY910861) was used as a MHV wild-type 
control. Delayed brain tumor cells (DBT-9), which are naturally permis- 
sive for MHV infection, and baby hamster kidney 21 cells expressing the 
MHV receptor (BHK-MHVR) were used for all experiments (48). Dul- 
becco’s modified Eagle medium (DMEM) (Gibco) was supplemented 
with 10% heat-inactivated fetal calf serum (FCS) for all experiments de- 
scribed. BHK-MHVR medium was supplemented with G418 (0.8 mg/ml; 
CellGro) to maintain selection for MHVR expression. All biochemical 
experiments were carried out using rabbit polyclonal antibodies previ- 
ously described in the literature. The antisera used include nsp2-specific 
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TABLE 1 Oligonucleotide primers used in the mutagenesis of MHV 
cDNA C fragment plasmids 


Mutation(s) made Primer sequence (5’ — 3’)* 


$133A ACG CTT CGT AGT GCC CAT ACC ATA AAG 

$133A/H134Y ACG CTT CGT AGT GCC TAT ACC ATA AAG 

$133A/T129M TTC CAT GTT ATG CTT ATG CTT CGT AGT 
GCC CAT ACC ATA 

$133A/H134Y/T129M TTC CAT GTT ATG CTT CGT AGT GCC TAT 
ACC ATA AAG 

H134A CTT CGT AGT AGC GCT ACC ATA AAG GGC 

H134Y ACG CTT CGT AGT AGC TAT ACC ATA AAG 
GGC TCC 

F219L AAC AGA TGC AAC TGG CTA GTG CAA AGT 
GAT AGT 

H270HH GCT ATT AAG AGG CTG CAT CAT TCT GGA 
TTC CAG 

E285V GGT AGT TGT GTG CTT GTT GAT GAG ACA 
CCA AGT 


* The nucleotide changes are underlined. 


(VU154), nsp5-specific (VU6), and nsp8-specific (VU123) antibodies (8, 
27, 39). 

Recovery of MHV mutant viruses. The MHV nsp5 mutant viruses 
were engineered through the infectious cDNA assembly strategy de- 
scribed previously by Yount et al. (48). In brief, the seven cDNA fragments 
were digested, gel purified, and ligated overnight at 16°C. Transcription of 
the extracted ligated DNA, as well as N cDNA which encodes the nucleo- 
capsid protein, was performed by using the mMachine T7 transcription 
kit (Ambion) under conditions previously described in detail (48). The 
transcribed genome and N gene were electroporated into BHK-MHVR 
cells, and the electroporated cells were placed into a subconfluent flask of 
DBT-9 cells and incubated at either 30°C for potentially ts viruses or 37°C. 

Mutagenesis of MHV cDNA C fragment. Assembly of the complete 
MHV genome is generated through the ligation of seven cDNA fragments 
(A to G) digested from individual plasmids as previously described by 
Yount et al. (48). All viruses were engineered by inserting the specified 
amino acid substitution into the MHV infectious clone (MHVic) C frag- 
ment containing the nsp5 sequence, which was constructed by PCR and 
cloned into the XL-pSMART vector (48). Sense and antisense primers 
were designed to be overlapping with nucleotide changes in the middle of 
the primer. The primers used for mutagenesis are listed in Table 1. The 
sequences of all mutant plasmids were confirmed prior to ligation and 
MHVic assembly. 

RNA extraction and sequencing. A confluent monolayer of DBT-9 
cells in a T25 flask was infected with viral mutant stocks at a multiplicity of 
infection (MOI) of 10 PFU per ml and grown until approximately 30 to 
50% of the cells had formed syncytia. Supernatant was removed from each 
T25 flask containing isolated mutant virus and stored at —20°C. The cells 
were harvested in TRIzol reagent (Invitrogen) for isolation of total RNA. 
Reverse transcriptase PCR (RT-PCR) was performed using SuperScript 
III RT (Invitrogen) and random hexamers (Applied Biosystems) at 55°C 
for 1 h, and the resulting cDNA was PCR amplified using Easy-A high- 
fidelity PCR cloning enzyme (Stratagene) and MHV genome oligonucle- 
otides covering the nsp5 region. Amplified regions were gel purified and 
analyzed by sequencing. 

Isolation and expansion of suppressor mutants. Confluent mono- 
layers of DBT-9 cells were initially infected with the temperature-sensitive 
viruses and incubated at either 40°C (S133A) or 30°C (F219L) with in- 
creases in temperature to 37°C and subsequently to 40°C. A plaque assay 
was then performed using 10-yl and 100-yl aliquots of viral stocks. 
DBT-9 cells were infected in duplicate using 6-well plates with serial dilu- 
tions of potential revertant samples with a 1-h adsorption period. The 
overlay contained a 1:1 mixture of 2% agar and 2X Dulbecco’s modified 
Eagle medium. The plates were incubated for 48 h so the plaques were 
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easily visible. Ten plaques were picked for each virus sample and resus- 
pended in gel saline. Isolated plaques for each virus were used to infect 
separate T25 flasks for expansion at 40°C. The flasks were removed from 
nonpermissive temperatures when 70 to 95% of cells had formed syncytia, 
and then RNA was isolated as described above. 

Analysis of viral growth and nsp5 processing. Confluent monolayers 
of DBT-9 cells in 60-mm dishes were infected at an MOI of 1 PFU/cell for 
growth analysis or 5 PFU/cell for immunoprecipitation analysis. In the 
temperature shift experiments, the cells were shifted from 30°C to 40°C at 
6 h postinfection (p.i.). During the growth analysis, samples of superna- 
tant were acquired and prewarmed medium was added back to maintain 
a fixed volume of medium for the cells. Virus titers were determined by 
plaque assay in duplicate. Immunoprecipitation experiments were carried 
out as previously described (40). Eluted proteins were resolved by sodium 
dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) on 4 to 
12% polyacrylamide gradient Bis-Tris gels (Nu-PAGE; Invitrogen) and 
analyzed by autoradiography. A ‘C-labeled high-molecular-weight stan- 
dard (New England BioLabs) and a full range rainbow marker (RPN 800E; 
Invitrogen) were used as protein mass standards. 

Sequence alignments and modeling of MHV nsp5 structures. The 
X-ray crystal structure of HCoV-HKUI nsp5 (Protein Data Bank [PDB] 
code 3D23) was used as a structural model of comparison (49). Structural 
models were generated using Modeler (16) and MacPyMol (DeLano Sci- 
entific). Other CoV nsp5 X-ray structures used for alignment and com- 
parison were SARS (PDB code 2H2Z), HCoV-229E (PDB code 1P9T), 
and IBV (PDB code 2Q6D) (2, 46, 47). A query of nsp5 sequences in 
GenBank resulted in the identification of 432 complete coronavirus nsp5 
amino acid sequences. A perl script was used to generate a representative 
sequence from repetitive sequences and eliminate redundancies as a 
means of unbiasing the data set; a final alignment of 130 unique, nonre- 
dundant sequences was prepared. A model of 100% conserved residues 
was determined using Consurf (4) and a sequence logo was generated 
using the WebLogo server (13) with the prepared nonredundant CoV 
nsp5 sequence data set. Distance calculations were determined by mea- 
suring from alpha carbon to alpha carbon in the structure of HCoV- 
HKUI nsp5. 


RESULTS 


Domain 2 $133A is a novel temperature-sensitive mutant of 
MHV nsp5. We previously identified two second-site mutations 
that suppressed the temperature-sensitive phenotype of MHV 
tsV148A: H134Y and S133N (40). Both H134 and S133 are greater 
than 20 A from V148, from the active site cavity, and from the S1 
substrate-binding site. The mechanism by which these residues 
complement or suppress tsV148A remains unclear. Further, sub- 
stitution of the H134Y or S133N residue alone in the isogenic 
cloned wild-type (WT) MHV background had little or no effect on 
virus replication at 37°C. We therefore sought to determine 
whether these residues had critical roles in nsp5 activity by engi- 
neering alanine substitution mutations in the MHV genome. 
S133A and H134A mutant viruses were recovered at 30°C and 
sequencing from the initial passage (passage 1 [P1]) 30°C virus 
stocks confirmed the presence of mutation $133A or H134A, with 
no other sequence changes in nsp5. The viruses were compared 
with recombinant WT MHV and with tsV148A for virus titer at 
30°C and 40°C, and the efficiency of plating (EOP) was calculated 
(titer at 40°C/titer at 30°C) (Fig. 2A). WIT MHV had an EOP of 3.3, 
a 3.3-fold increase in visible plaques at 40°C compared to 30°C. 
The previously described tsV148A virus exhibited an EOP of 3 X 
10~°, confirming the ts phenotype reported by Sparks et al. (40). 
The engineered H134A mutant had an EOP of 101, which is 
similar to that of the previously reported engineered H134Y mu- 
tation in the WT MHV background (40). In contrast, the $133A 
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FIG 2 Analysis of replication of tsS133A and suppressor mutants. (A) Virus 
titers and EOP data for WT MHV and S133A mutant viruses determined by 
plaque assay at 30°C and 40°C. The virus titers were calculated in duplicate by 
infection of DBT-9 cells. Recombinant second-site suppressor mutants (*) 
and engineered recombinants whose mutations were artificially recombined 
(A) are indicated. (B and C) Growth analysis of the WT MHV, recombinant 
tsS133A, suppressor mutants (S133A/H134Y and $133A/T129M/H134Y), and 
engineered mutant (S133A/T129M) grown at 30°C (B) or grown at 30°C and 
then shifted to 40°C at 6 h p.i. on DBT-9 cells (C). The virus titers were 
determined by plaque assay on DBT-9 cells at 30°C. 


mutant had an EOP of 2.7 X 10°, consistent with an independent 
ts phenotype (Fig. 2A). Interestingly, reversion of tsV148A had 
resulted in identification of an allele at $133 that with an Asn 
substitution could suppress the V148A ts phenotype (S133N), but 
with a conservative Ala substitution resulted in a new ts phenotype 
(S133A). H134A was not further examined, and tsS133A was stud- 
ied in subsequent experiments. 

Phenotypic reversion of tsS133A at the nonpermissive tem- 
perature reveals second-site suppressor mutations in domain 2. 
To test for revertants or suppressors of S133A, DBT-9 cells were 
infected with tsS133A at 37°C, followed by a shift to 40°C. Recov- 
ered supernatant virus was used to infect cells at 40°C, and 10 virus 
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plaques were isolated and expanded at 40°C. RT-PCR sequencing 
of the nsp5 coding regions from all 10 clones confirmed the reten- 
tion of the $133A (j9¢9;AGC to GCC with the mutation shown 
underlined) engineered mutation, as well as two distinct patterns 
of second-site nonsynonymous mutations resulting in amino acid 
substitutions: H134Y (,o69gCAT to jo69g AT) and T129M/H134Y 
(,0503ACG to 19593 A LG; jg60gCAT to i960 LAT). No primary rever- 
sion at S133A was identified, and no other mutations were de- 
tected in the nsp5 coding sequence. The $133A/H134Y plaque 
isolate had an EOP of 0.5, and the $133A/T129M/H134Y plaque 
isolate had an EOP of 2, suggesting almost complete suppression 
of the S133A ts phenotype (Fig. 2A). To test the contribution of the 
second-site mutations to suppression of the ts phenotype, we en- 
gineered different combinations of mutations into the isogenic 
cloned MHV background. Since both biologically derived sup- 
pressor mutants contained H134Y and there were subtle differ- 
ences in the EOP values, we tested the independent contribution 
of the T129M substitution by introducing the $133A and T129M 
mutations in the absence of the H134Y mutation. All of the engi- 
neered recombinant mutants were readily recovered at 30°C, and 
sequencing confirmed that the engineered changes were present 
and no other mutations had arisen in nsp5. The EOPs of the re- 
combinant $133A/H134Y and $133A/T129M/H134Y mutants 
were identical to the cognate biologically recovered mutant, dem- 
onstrating that the identified changes in nsp5 were necessary and 
sufficient for the phenotypic reversion (Fig. 2A). In contrast, the 
engineered recombinant $133A/T129M mutant showed an EOP 
of 10%, 2 log units greater than ts$133A but still significantly ts 
compared to either the $133A/H134Y or S133A/T129M/H134Y 
mutant. The results demonstrated that H134Y was sufficient for 
suppression of the tsS133A, T129M was unable to suppress the ts 
phenotype by itself, and the combination of T129M and H134Y 
was additive or synergistic, suggesting that they are two distinct 
mutations that could have arisen sequentially or concurrently in 
the same or different genomes. 

Growth of tsS133A and second-site suppressor mutant vi- 
ruses demonstrate differential effects of single and multiple 
suppressor mutations. The capacity to generate a visible plaque is 
one indicator of viral fitness but has limited ability to predict over- 
all growth fitness. We therefore compared the capability of WT 
and mutant viruses to replicate at permissive (30°C) and nonper- 
missive temperatures (40°C) in single-cycle growth experiments. 
On the basis of earlier single-cycle growth studies of MHV repli- 
cation, we infected replicate plates of DBT-9 cells with WT and 
recombinant mutant viruses at 30°C and a multiplicity of infec- 
tion (MOI) of 1 PFU/cell (18, 19). At6h p.i., one replicate plate for 
each virus was transferred to the nonpermissive temperature of 
40°C and one plate was maintained at 30°C, with supernatant 
samples obtained from 0 to 30 h p.i. at regular intervals for deter- 
mination of viral titers. WT virus incubated at 30°C demonstrated 
onset of exponential growth between 10 and 12h p.i. and achieved 
peak virus titers of ~ 10° PFU/cell by 30 h p.i. (Fig. 2B), consistent 
with previous studies (40). In contrast, tsS133A and all three re- 
combinant suppressor mutant viruses demonstrated identical 
growth curves, with a 1-log-unit reduction in titer compared to 
the titer of WT MHV from 12 to 24 h p.i., but achieving titers 
identical to WT virus by 30 h p.i. 

When the replicate virus-infected monolayers were shifted 
from 30°C to 40°C at 6 h p.i. (Fig. 2C), the WT virus-infected 
monolayers showed onset of exponential growth within 2 h after 
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shift and achieved peak titers between 12 and 16h p.i, albeit at 
lower titers due to rapid destruction of the monolayer. The 
tsS133A mutant virus showed profoundly impaired growth for 10 
h after the temperature shift. Both recombinant $133A/H134Y 
and §133A/T129M/H134Y mutant viruses demonstrated onset of 
exponential growth and peak titers similar to those of the WT 
virus following the shift, consistent with EOP analysis. However, 
the S133A/T129M virus showed a significant delay in exponential 
growth before achieving titers similar to those of the WT virus by 
6 h after the temperature shift. These results collectively indicate 
that the S133A mutation confers a subtle replication defect at 30°C 
that is not further impaired nor complemented by the suppressor 
mutations and that H134Y alone or in combination with T129M is 
sufficient for suppression of the tsS133A growth phenotype. The 
S133A/T129M virus, while demonstrating improved growth over 
tsS133A, was still impaired compared to mutants containing 
H134Y, a result consistent with the EOP data, and supporting 
either conjecture that T129M arose first to be superseded by emer- 
gence of H134Y or arose second by conferring a subtle growth 
advantage of the combination over H134Y alone. 

Recovery and reversion of recombinant nsp5 tsF219L. The 
experiments with tsS133A, in combination with our previous 
studies, demonstrated that the H134Y substitution was able to 
suppress two distinct and independently derived nsp5 ts alleles, 
tsV148A and tsS133A, suggesting an important role for intra- or 
intermolecular communication involved in regulation of nsp5 ac- 
tivity. However, both of these ts alleles are in domain 2 as were the 
suppressor mutations, and thus could not provide insight into 
potential long-distance communication between domains. We 
therefore next sought to determine whether H134 would emerge 
as a suppressor for another putative ts F219L allele in domain 3. 
The 9364 UUU-to-CUU mutation resulting in the F219L substitu- 
tion was predicted as a ts allele by Sawicki et al. using partial ge- 
nome sequencing and reversion analysis of the ts mutant Alb ts16 
(36). However, this was not confirmed in this study as the sole 
mutation by complete genome sequencing or by reverse genetic 
analysis. In addition, primary reversion occurs rapidly in biolog- 
ical mutants with a single-nucleotide polymorphism; therefore, 
possible second-site suppressors could not be identified. We en- 
gineered the F219L codon change as a two-nucleotide mutation 
(,0864UUU to CUA) in the isogenic MHV clone, which would 
require a two-nucleotide change for reversion to Phe219 (UUU or 
UUC). The engineered recombinant F219L mutant was recovered 
at 30°C, and complete genome sequencing confirmed the 
10864 UUU-to-CUA mutations as the only changes in the genome. 
The recombinant F219L mutant had an EOP of 3 X 107°, con- 
firming that the F219L substitution alone was sufficient to confer 
a ts phenotype (Fig. 3A). The titer and plaque morphology of the 
tsF219L mutant were indistinguishable from those of WT MHV at 
30°C (data not shown). 

Identification of second-site suppressor mutations of recom- 
binant tsF219L. To select for phenotypic revertants or suppressor 
mutations of tsF219L, DBT-9 cells were infected with recombi- 
nant fsF219L mutant virus at 40°C. However, no cytopathic effect 
(CPE) or productive infection occurred at 40°C, despite multiple 
attempts and prolonged incubation. Consequently, we initiated 
infection at 30°C for 6 h, followed by a shift to 37°C for 24 h. This 
stock was then passaged at 37°C with a shift to 40°C, followed by 
passage and selection of 10 plaques at 40°C. Sequencing of 10 
plaque clones confirmed retention of the engineered jog¢4CUA 
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FIG 3 Analysis of replication of tsF219L and second-site suppressor mutants. 
(A) Virus titers and EOP data for WT MHV and F219L mutant viruses deter- 
mined by plaque assay at 30°C and 40°C. Titers were calculated in duplicate by 
infection of DBT-9 cells. Recombinant second-site suppressor mutants (*) 
and engineered recombinants whose mutations were artificially recombined 
(A) are indicated. (B and C) Growth analysis of WT MHV, recombinant 
tsF219L, suppressor mutants (F219L/H134Y/E285V and F219L/H134Y/ 
H270HH), and engineered mutants (F219L/H270HH/E285V, F219L/E285V, 
F219L/H270HH, and F219L/H134Y) grown at 30°C (B) or grown at 30°C and 
then shifted to 40°C at 6 h p.i. (C) on DBT-9 cells. The virus titers were 
determined by plaque assay on DBT-9 cells at 30°C. 


(F219L) codon. However, all 10 isolated plaque cloned viruses dem- 
onstrated one of two patterns of second-site mutations in nsp5 in the 
presence of tsF219L: (i) H134Y (;969gCAT to j9¢60gLAT)/H270 dupli- 
cation (H270HH) (1491:6CAT to 1 49;6CATCAT) with an EOP of 1.2 or 
(ii) H134Y (CAT to TAT)/E285V (1196GAA to 449¢,GTT) with an 
EOP of 2.9. To test the contribution of the H134Y, E285V, and 
H270HH changes to suppression of tsF219L, the identified mutations 
were engineered with and without F219L, alone or in the combina- 
tion seen in the recovered viruses. In addition, although no F219L/ 
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E285V/H270HH mutant was identified among the sequenced 
plaques, we engineered this combination to test for the capacity of the 
combination to also suppress the ts phenotype. In total, nine genomes 
were engineered: F219L/E285V, F219L/H134Y, F219L/H270HH, 
F219L/H134Y/E285V, F219L/H134Y/H270HH, F219L/H270HH/ 
E285V, H134Y/E285V, H134Y/H270HH, and H270HH/E285V. All 
viruses were recovered at 30°C with the engineered mutations 
detected and confirmed by sequencing across nsp5. The EOP was 
determined for each of the viruses by plaque assay on DBT-9 cells 
at 30°C and 40°C (Fig. 3A). The second-site substitution combi- 
nations identified by reversion analysis were sufficient to suppress 
the F219L ts phenotype: H134Y/E285V, EOP = 2.9; and H134Y/ 
H270HH; EOP = 1.2. When the second-site suppressor alleles 
were tested alone with F219L, the results showed that the individ- 
ual substitutions either minimally or partially suppressed tsF219L: 
H134Y, EOP = 2 X 10°; E285V, EOP = 1 X 10 °;and H270HH, 
EOP = 2 X 10 *. When the second-site substitutions were intro- 
duced in the WT background, either alone or in combination, 
there was no effect on EOP, suggesting that the changes were not 
responsible for any replication defects in the presence of F219L. 
The nonbiologically derived combination of F219L/E285V/ 
H270HH also completely suppressed the ftsF219L phenotype 
(F219L/E285V/H270HH, EOP = 1.0). The EOP results confirmed 
that the biologically identified second-site substitution combina- 
tions were both necessary and sufficient to suppress the F219L ts 
phenotype. Further, the results showed that H134Y did emerge as 
a suppressor allele for domain 3 tsF219L, but that in contrast to 
tsV 148A and tsS133A, suppression of tsF219L required at least one 
other substitution in domain 3 in combination with H134Y. 

Growth of tsF219L and revertants. The tsF219L and recombi- 
nant suppressor mutant viruses were grown in DBT-9 cells at 30°C 
or beginning at 30°C with a shift to 40°C at 6 h p.i. (Fig. 3B and C). 
At 30°C, the growth kinetics and virus yields for all mutant viruses 
were indistinguishable from those of WT MHV. Following the 
shift to 40°C at 6 h p.i., fsF219L showed no further replication for 
10 h. In contrast, the mutant viruses containing any two of the 
H134Y, E285V, and H270HH substitutions showed growth after 
the shift to 40°C similar to WT. In contrast, all single suppressor 
residues expressed with F219L demonstrated a 4-h lag before ex- 
ponential growth compared to the double mutants, but ulti- 
mately, they achieved peak titers similar to those of the double 
mutants. The titers at 10 h p.i. were consistent with the EOP data 
(10~* to 10-* compared to WT) and overall demonstrated that 
the individual mutations were capable of improved viral replica- 
tion compared to the tsF219L virus yet still were ts compared to 
WT (Fig. 3A and C). Thus, in contrast to tsV148A and tsS133A, 
suppression of tsF219L appears to require a combination of at 
least two second-site mutations. This may be the explanation for 
the tight tsF219L phenotype at 40°C, as well as for the necessity of 
sequential passage of the tsF219 mutant virus at 30°C, 37°C, and 
40°C to recover phenotypic revertants. In addition, these results 
identify H134Y as a second-site suppressor for a third fs allele in 
nsp5, this one in a domain 3. Overall, the results demonstrate 
cooperation of H134, E285, and H270 in nsp5 for efficient virus 
replication, as well confirming communication between nsp5 do- 
mains 2 and 3. 

The tsS133A and tsF219L mutant viruses have temperature- 
sensitive impaired processing by nsp5. The nsp5 protease is pres- 
ent in all CoVs and is responsible for 11 maturation cleavage 
events in the replicase polyprotein (nsp4 through nsp16). To di- 
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rectly compare the nsp5 protease activity of tsV148A, tsS133A, 
tsF219L, and second-site suppressor mutants, DBT-9 cells were 
infected with WT and mutant viruses at an MOI of 5 PFU/cell and 
incubated at 30°C. At 6h p.i., replicate monolayers were main- 
tained at 30°C or transferred to 40°C, and infected cells were ra- 
diolabeled with [°°S]Met-Cys. Lysates from infected, radiolabeled 
cells were immunoprecipitated with antibodies specific for nsp2, 
nsp5, and nsp8 to test for processing of nsp2 by PLP1 and of nsp5 
and nsp8 by nsp5. 

Immunoprecipitation with antibodies specific for nsp2 of all 
lysates from WT and mutant virus-infected cells with labeling at 
30°C and 40°C resulted in detection of mature processed nsp2, 
demonstrating that at both permissive and nonpermissive tem- 
peratures, there was similar translation of ppla (nsp1 to nsp11) 
and normal PLP1 activity (Fig. 4). Recently, Stokes et al. reported 
that a ts mutation in nsp3 resulted in a significant decrease in 
nsp5-mediated processing (41). In our study, we detected the 
nsp4-to-nsp10 precursor polyprotein (150 kDa), which demon- 
strates that PLP2 is functional and is not inhibited by the nsp5 
mutations (Fig. 4). The presence of mature nsp2 and the p150 
bands at the nonpermissive temperature indicate that both PLP1 
and PLP2 domains of nsp3 are active and are not affected by the 
nsp5 mutations. Although we did not test the processing of nsp3 
directly, detection of both p150 and mature nsp2 is consistent 
with normal processing of N and C termini of nsp3. Immunopre- 
cipitation of cells infected with WT MHV by nsp5-specific anti- 
bodies at both 30°C and 40°C resulted in detection of mature 
processed nsp5, as well as coimmunoprecipitation of nsp8, and 
two distinct small protein bands at 10 and 12 kDa, which is con- 
sistent with the predicted migration of nsp7 and nsp9. Immuno- 
precipitation with nsp8-specific antibodies detected nsp8 as well 
as probable coprecipitation of nsp5 and the 10- and 12-kDa pro- 
teins. These results show that expression and processing of ppla 
proteins nsp5, nsp7, nsp8, and nsp9 are accelerated at 40°C in cells 
infected with WT MHV. 

The temperature-sensitive viruses tsS133A and tsF219L, as well 
as the previously described tsV148A, exhibited profoundly im- 
paired processing of nsp5 and nsp8 at 40°C compared to 30°C, 
indicating a specific defect in processing by the nsp5 protease. This 
was consistent with decreased detection of the reciprocal coim- 
munoprecipitating protein (nsp8 or nsp5) as well as decreased 
detection of the 10- and 12-kDa proteins. Although the detection 
of nsp5-processed proteins was profoundly decreased, we were 
unable from multiple replicate experiments (more than 5 experi- 
ments) to demonstrate complete loss of nsp5 activity. The results 
suggest that the $133A, F219L, and V148A mutations do not di- 
rectly affect the catalytic or substrate-binding functions of nsp5 
but rather modify protease activity in other ways. Alternatively, 
it is possible that the residual processing might be the result of 
nsp5 expressed and folded into active forms or complexes prior 
to the temperature shift, and thus, the protein still retains re- 
sidual activity. 

Viruses carrying biological and engineered suppressor muta- 
tions demonstrated restoration of processing by nsp5 that directly 
correlated with the degree of recovery of EOP and virus growth. 
Single second-site mutants $133A/T129M, F219L/H134Y, F219L/ 
H270HH, and F219L/E285V showed an increase in detectable 
processed nsp5 and nsp8 only, while the double second-site sup- 
pressors restored WT-like patterns of processed proteins. Collec- 
tively, the results show a direct correlation of detection of proteins 
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FIG 4 Proteolytic processing of WT, tsS133A, tsF219L, and suppressor viruses. DBT-9 cells were *°S radiolabeled during viral infection or mock infection. 
Cellular lysates were harvested from cells infected with WT MHV, tsS133A mutants, tsF219L mutants, and previously described tsV 148A and from mock-infected 
cells (40). Labeled proteins were immunoprecipitated using antiserum specific for nsp2, nsp5, and nsp8. The temperature during virus infection is indicated 
above the lanes (30°C [30] or a temperature shift from 30° to 40°C at 6 h p.i. [40]). One hundred microliters of lysate was used for all immunoprecipitations. The 
positions of putative viral proteins are shown to the right of the gels based upon the predicted size, and identified bands are labeled. The positions of molecular 
weight standards (MW) are shown to the left of the gels, and sizes are shown in kilodaltons. Protein expression profiles were resolved by SDS-PAGE and visualized 


by autoradiography. a-nsp2, anti-nsp2 antibody. 


processed by nsp5 (nsp5 and nsp8) and the extent of restored 
growth fitness in culture and EOP analysis in plaque assay. Fur- 
ther, these findings indicate that the impairment in growth at 
nonpermissive temperatures is not due to complete inactivation 
of nsp5 protease activity. 

Analysis of ts and suppressor mutations in nsp5 structures. 
To evaluate the structural relationships between the ts and sec- 
ond-site suppressor mutations, the distances between the combi- 
nations of ts and suppressor residues were determined by analyz- 
ing the crystal structure of the nsp5 protease of a closely related 
human CoV, HKU1. The structure of MHV nsp5 has not been 
determined; however, MHV and HKU1 nsp5 proteases exhibit 
84% sequence identity (Fig. 1C) and share all of the same amino 
acids at the residue positions reported in this study with the ex- 
ceptions of H134 and H270 (Y134 and Y270 in HKU1). All of the 
second-site suppressor residue positions in HKU1 nsp5 were 
greater than 10 A distance from the ts residues in the monomeric 
structure with the sole exception that the H134 residue is 3.8 A 
from the juxtaposed S133 (data not shown). Measurement of the 
distance between residues in different monomers of the dimeric 
structure of both SARS-CoV and HKU1 nsp5 demonstrated that 
no two residues from this study were closer than 15.0 A. Measure- 
ment of the distance between the ts and second-site suppressor 
mutants and the catalytic dyad residues, H41 and C145, showed 
that only V148A was within 10 A of either residue. Collectively, 
these data demonstrate that the relationship between the ts and 
suppressor mutations cannot be explained by direct interactions 
between residues and that nsp5 dimerization does not provide 
direct intermonomer associations between the residues identified 
in this study. Modeling of the $133A and F219L mutations on the 
structure of HKU1 nsp5 failed to predict any clear pathways of 
side chain remodeling or perturbation between the ts residues and 
the protease active site (data not shown). In contrast, analysis of 
residue conservation using an alignment of 130 nonredundant 
CoV nsp5 amino acid sequences identified a series of 100% iden- 
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tical residues that span the regions of nsp5 between each of the ts 
residues (S133, V148, and F219) and the common second-site 
suppressor residue (H134) (Fig. 5). These findings indicate that 
the structural and functional perturbations on nsp5 protease of 
the ts mutations may span long distances across the protease struc- 
ture through yet to be identified cooperative interactions. 
Residue conservation of ts and suppressor alleles is group 
specific. To evaluate the variability at the ts and suppressor alleles, 
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FIG 5 Conserved, ts, and suppressor alleles in a solved CoV nsp5 structure. A 
conservation map of 100% identical nsp5 residues (black) across 130 nonre- 
dundant CoV nsp5 sequences is shown on the HKU1 nsp5 protease monomer 
structure. The locations of identified ts alleles (red) and suppressor alleles 
(green) are indicated. 
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FIG 6 Coronavirus protein sequence conservation of ts and suppressor alleles. 
(A) A sequence logo of conservation of ts and suppressor mutations across an 
alignment of 130 nonredundant CoV nsp5 sequences was generated using 
WebLogo (13). The height of each letter corresponds to the relative conserva- 
tion of that amino acid at the position, and the height of the column corre- 
sponds to the sequence conservation at the position. The residue numbers are 
relative to the MHV amino acid positions. (B) The ts and suppressor residues 
of coronavirus species by CoV genus are shown for each MHV residue posi- 
tion. a-CoV1, alphacoronavirus 1; TGEV, transmissible gastroenteritis virus; 
PEDV, porcine epidemic diarrhea virus; BtCoV, bat coronavirus. 


a sequence logo (13) was generated using the nonredundant align- 
ment of 130 CoV nsp5 amino acid sequences, and the ts and sup- 
pressor alleles for 17 coronavirus species were analyzed (Fig. 6A 
and B). From two to eight different residues occupy each position 
across available CoV sequences (Fig. 6A). Surprisingly, the H134Y 
common second-site suppressor mutation selected for a tyrosine 
that is already present in several betacoronaviruses. Further, the 
S133N second-site suppressor first reported by Sparks et al. is 
common as an Asn in many coronaviruses (40). Conservation of 
distinct ts and suppressor alleles within the three genera suggest 
that there may be select combinations of alleles that are necessary 
for nsp5 activity (Fig. 6B). These findings further suggest that 
alterations in structure attributed to these residues could have 
analogous combinations in other CoVs. 


DISCUSSION 


Previous studies with other viruses as diverse as HIV, Sindbis vi- 
rus, poliovirus, and vaccinia virus have reported ts mutations in 
virus proteases that affect protein processing, RNA synthesis, and 
virus capsid assembly (3, 23, 24, 29). Most of the ts alleles of de- 
scribed viral proteases have occurred at individual conserved res- 
idues or in pairs of conserved and structurally adjacent residues. 
For the HIV protease, residues distant from known catalytic and 
functional determinants have been shown to be critical for pro- 
tease activity (22, 34). The current study of MHV nsp5 (3CLpro) 
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extends our understanding of RNA virus proteases by demon- 
strating that multiple nonconserved and structurally distant resi- 
dues in the CoV nsp5 protease participate in long-distance com- 
munication within and between the protease structural domains, 
function cooperatively to suppress ts phenotypes, and are impor- 
tant for nsp5 activity during replicase polyprotein processing. 
Specifically, we demonstrate that independent and physically dis- 
tant ts alleles in domains 2 and 3 resulted in selection of the same 
H134Y ts suppressor allele, as well as selection of several additional 
suppressor alleles. With the exception of tsV148A, all of the resi- 
dues are structurally distant from the catalytic, substrate, and 
dimerization residues in monomers and between nsp5 molecules 
in solved nsp5 dimer structures. In addition, suppressor muta- 
tions arose during reversion analysis in combinations that were 
required for or augmented restoration of nsp5 processing. Fur- 
ther, artificial combinations of suppressor mutations not seen 
during biological reversion analysis also showed cooperative sup- 
pression of the ts phenotype. These results all support the hypoth- 
esis that long-distance communication occurs between multiple 
residue nodes to regulate nsp5 activity. 

Since nsp5 must recognize and process 11 closely related cleav- 
age sites in the setting of rapidly changing substrate and cleavage 
site concentrations from polyproteins and processing intermedi- 
ates, any explanation for nsp5 intra- or intermolecular communi- 
cation needs to account for these evolving variables. The observa- 
tion that a series of residues that are completely conserved across 
the coronaviruses span the regions between the ts and second-site 
mutations in solved and modeled CoV nsp5 structures (Fig. 5) 
suggests a possible mechanism for such communication. The po- 
tential linkage by conserved residues in the nsp5 tertiary structure 
would be similar to the model proposed by Ranganathan and co- 
workers in which networks of distant residues are connected in 
tertiary structure and are nodes for allosteric communication (22, 
42). These networks may further demonstrate coevolution of res- 
idues that maintain the protease structure. 

Temperature-sensitive mutations and nsp5 activity. Analy- 
ses of the MHV nsp5 ts mutant viruses demonstrated a defect in 
nsp5-mediated processing and virus replication at the nonpermis- 
sive temperature. However, mature nsp5 cleavage products were 
still detected for all three ts mutant viruses at the nonpermissive 
temperature, albeit at profoundly reduced levels, indicating that 
defects in viral growth at the nonpermissive temperature were not 
the result of complete loss of nsp5 activity. We previously reported 
that nsp5 activity was not present in tsV148A at the nonpermissive 
temperature, but in that study we did not test nsp5 processing at 
cleavage sites flanking nsp5. In the present study, we determined 
that tsV148A functions like tsS133A and tsF219L, exhibiting more 
impaired processing of nsp8 than nsp5 at the nonpermissive tem- 
perature, while retaining residual nsp5 activity at the flanking 
nsp4-nsp5 and nsp5-nsp6 cleavage sites (40). Overall, our results 
suggest that the nsp5 ts phenotype may be due to an alteration of 
protease activity at different polyprotein substrate sites. Altered 
cleavage site specificity has been reported for mutations intro- 
duced in the interdomain loop between domains 2 and 3 of the 
arterivirus equine arteritis virus (EAV) nsp4 protease, a structural 
orthologue of CoV nsp5 (43). Several of the recovered EAV nsp4 
interdomain loop mutants did not abolish nsp4 protease activity 
but rather altered the substrate specificity. It is possible that the ts 
mutants of nsp5 are acting in a similar manner. 

Nonconserved residues in nsp5 protease function. Align- 
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ment of 130 nonredundant CoV nsp5 sequences demonstrates 
that the residues that resulted in the nsp5 ts mutations are not 
conserved. There was more conservation of specific residues 
within the betacoronavirus genus, including MHV and SARS- 
CoV, specifically at residues V148, F219L, and E285. Finally, the 
revertants H134Y and $133N selected amino acids that already 
were present in nsp5 sequences of other alpha- and betacoronavi- 
ruses. Among the second-site residues, the H270 position shows 
the greatest variation with eight different amino acid residues tol- 
erated across the CoVs. The tolerance for a wide variety of residues 
at this position may explain the acceptance of a duplicated codon 
resulting in a second histidine residue at the position. While the 
results have to be interpreted in light of the overall nsp5 sequence 
variability or conservation, they suggest that the ts and revertant 
residues may represent an evolutionarily adaptive network co- 
evolving with other interacting proteins or cleavage sites. In sup- 
port of this hypothesis, it was reported that a mutation at the P1 
position in the cleavage site between nsp15 and nsp16 of the CoV 
infectious bronchitis virus resulted in a debilitated virus whose 
phenotype was compensated for by a mutation in nsp5 (17). Sim- 
ilarly, ts mutations in MHV nsp3 and nsp10 resulted in altered 
nsp5-mediated nsp5 processing (14, 41). A network of noncon- 
served and mutationally flexible residues could account for the 
rapid emergence of second-site revertants of ts viruses with com- 
binations of more than one second-site mutation that also func- 
tion in combinations not derived during virus reversion analysis. 

Models for testing nsp5 long-distance communication. 
Coronavirus nsp5 functions in the setting of the largest known 
RNA virus polyprotein and must orchestrate 11 distinct cleavage 
events. It is clear from our findings and others that nsp5 activity is 
affected by changes within the structure of the protease, by 
changes at nsp5 cleavage sites, and by changes in other replicase 
proteins. Biochemical studies of the fs and revertant nsp5 mole- 
cules from this study should determine whether protein altera- 
tions and communication are intramolecular or within the nsp5 
dimer. Continued mutagenesis of nsp5 in the context of a virus 
will permit testing for allosteric interactions across the replicase 
polyprotein that impact specific ts and revertant alleles during 
virus replication. To test for additional communication nodes in 
nsp5, we will continue to use an iterative approach for mutagen- 
esis at ts and revertant alleles with alanine, nonconservative 
changes, and substitution of residues from other CoV nsp5 se- 
quences. Finally, since some of the MHV revertant alleles such as 
H134Y are already present in other CoVs, it may be possible to test 
whether the variability at the revertant residues across different 
coronaviruses affects whether a particular substitution results in a 
ts phenotype in different nsp5 backgrounds. The application of 
this data set in combination with structure information and bioin- 
formatic analysis should allow dissection of the extent and mech- 
anism of the communication network that regulates nsp5 activity 
during coronavirus replication. 
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